Zero-field Time Correlation Functions of Four Classical Heisenberg Spins on a Ring 



O 

o 

(N 



I 

C 

O 

o 



> 
o 
in 
o 
in 
o 



03 



i 

-a 
c 

o 
o 



Richard A. Klemm 1 ' 2 ' 3 and Marshall Luban 2 
1 Max-Planck-Institut fur Physik komplexer Systerne, Nothnitzer Strafie 38, D-01187 Dresden, Germany 
2 Ames Laboratory and Department of Physics and Astronomy, Iowa State University, Ames, IA 50011 USA 
3 Materials Science Division, Argonne National Laboratory, Argonne, IL 60439 USA 

(February 1, 2008) 



A model relevant for the study of certain molecular mag- 
nets is the ring of N = 4 classical spins with equal near- 
neighbor isotropic Heisenberg exchange interactions. Assum- 
ing classical Heisenberg spin dynamics, we solve explicitly for 
the time evolution of each of the spins. Exact triple inte- 
gral representations are derived for the auto, near-neighbor, 
and next-nearest-neighbor time correlation functions for any 
temperature. At infinite temperature, the correlation func- 
tions are reduced to quadrature. We then evaluate the Fourier 
transforms of these functions in closed form, which are double 
integrals. At low temperatures, the Fourier transform func- 
tions explicitly demonstrate the presence of magnons. Our ex- 
act results for the infinite temperature correlation functions in 
the long-time asymptotic limit differ qualitatively from those 
obtained assuming diffusive spin dynamics. Whether such ex- 
plicitly non-hydrodynamic behavior would be maintained for 
large- N rings is discussed. 

75.10.Hk,75.75.+a,75.30.Ds,75.75.-y 



I. INTRODUCTION 

Recently, there has been a rapidly growing interest 
in the physics of molecular magnets. [|l]j|] These com- 
pounds can be synthesized as single crystals of identi- 
cal molecular units, each containing several paramagnetic 
ions that mutually interact via Heisenberg exchange. The 
intermolecular (dipole-dipole) magnetic interactions are 
in the great majority of cases utterly negligible as com- 
pared to intramolecular magnetic interactions. Measure- 
ments of the magnetic properties therefore reflect those 
of the common, individual molecular units of nanometer 
size. Their dynamics can be studied by inelastic neutron 
scattering, as well as by nuclear magnetic resonance and 
electron paramagnetic resonance experiments. Some of 
these molecular magnets are made of very small clusters 
of magnetic ions. The smallest clusters are dimers of V 4+ 
(S = 1/2) and of Fc 3 + (S = 5/2), §|§ a nearly equilat- 
eral triangle array of V 4+ spins, H, a nearly square array 
of Nd 3+ (total spin j = 9/2), a regular tetrahedron of 
Cr 3+ (S = 3/2), JtJ a frustrated tetrahedral pyrochlore 
of Tb 3 + (5 = 5/2), g and a "squashed" tetrahedron of 
Fe 3+ spins. There has also been an example of a 

four-spin ring which is coupled to nearby rings, although 
the spin value (S — 1/2) is small, and thus requires a 



quantum treatment. JTT[ ] In addition, larger rings, most 
notably with 6, 8, or 10 Fe 3+ spins, have been studied. 

In some of these systems, the spin value of an indi- 
vidual magnetic ion is large enough that the dynamics 
can be closely approximated by the classical theory, as 
long as one does not go to temperatures that are too 
low. Thus, it is useful to study such systems theoreti- 
cally, in order to investigate the types of dynamical spin 
behavior that can occur. Such investigations can pro- 
vide helpful physical insight, as well as some guidance 
for systems that might be studied experimentally. It will 
also be interesting to compare the classical results with 
those emerging from studies of their quantum analogues, 
such as has been done for the dimer and the equilateral 
triangle. 

Perhaps more interesting, however, is the question as 
to whether the long-time asymptotic behavior of the two- 
spin correlation functions at infinite temperature will be 
consistent with the results of a hydrodynamic-like theory, 
in which the exact equations governing the spin dynam- 
ics are approximated by linear diffusion-like equations 
[see Eq. (jpj)]. This question has been the subject of 
much debate in the literature, [^7 18 H]j2(|, and a solu- 
tion of the spin dynamics for the four-spin ring might aid 
in our understanding of this more fundamental problem. 
Here we derive exact results which explicitly demonstrate 
that the infinite temperature, long-time asymptotic lim- 
its of the N = 4 two-spin correlation functions are non- 
hydrodynamic. 

The layout of the paper is as follows. In Sec. II, we 
give the notation, partition function, and derive the ex- 
act time evolution of the individual spin vectors. In Sec. 
Ill, we give the results for the time correlation functions. 
At infinite temperature, these results can be expressed 
as single integrals, but at finite temperatures, they are 
triple integrals. We also present our derivation of the 
Fourier transforms of the deviations of the correlation 
functions from their infinite time asymptotic limits. Fi- 
nally, we invite the reader to read our discussion and 
conclusions in Sec. IV, even if one has only a minimal 
interest in the mathematical developments presented in 
Sees. II and III. In this final section, we also discuss the 
non-hydrodynamic aspects of our exact results for the 
infinite temperature, long-time asymptotic behaviors of 
the two-spin correlation functions, and raise the question 
as to whether such non-hydrodynamic features might be 
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maintained for larger rings. 

II. SPIN DYNAMICS 

A. Notation and partition function 

We study the dynamics of four interacting spins on a 
ring. Each spin has unit magnitude and can assume an 
arbitrary direction. We suppose that it interacts only 
with its two nearest neighbors. We label these spins 
Si,S2, etc., where |Sj| = 1, illustrated in Fig. 1. 
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and a — (3J/2. Eq. @ was obtained previously, |||] 
and an analysis of the integral was also presented. The 
overall factor of ^ can be dropped, as it will not appear in 
the correlation functions. In deriving Z, we note that the 
angles S24 makes with k and q are independent, as are 
the angles S13 makes with k and p. Hence, it is helpful 
to first integrate over the angles Si and S3 make with p 
and the angles S2 and S4 make with k. Then, one may 
choose the coordinates of k and p relative to S24 and 
S13, respectively. 



FIG. 1. Sketch of a ring with four classical spins at the cor- 
ners, each interacting with its nearest neighbors with strength 
J. 

The Hamiltonian for this system is thus 
4 

H =-jJ2 s i- s i+i, (!) 

i=l 

= -(J/2)(S 2 -S 2 13 -S 2 24 ), (2) 
Si + S 3 , S 2 4 = S 2 + S4, and 



where S5 = Si, S13 
S = S 13 + s 24 . 

The partition function Z can be readily found from 
Eq. (0), leading to 



Z = Tr cxp(-/3iJ) 

flit 
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B. Exact time evolution 

The dynamics of the spins arise from the Heisenberg 
equations of motion, 



dSj 
~dt 



1 N 



x Sh 



(7) 



3=1 

<ij> 



where Sj^ + i = S; for any integer i. Our primary concern 
in this paper is the case N = 4, for which Eq. (@) may 
be rewritten as 



dS±3 

dt 
dS 2A 

dt 



-S1.3 x S 2 
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which lead to 



and 



dS 13 

dt 
dS 24 

dt 



dS 
~dt 



-S13 x S, 

T 

-S 24 x S, 

T 



0. 



(8) 
(9) 

(10) 
(11) 

(12) 



The phenomenological classical spin precession rate - 
can be obtained from first principles, starting from a 
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quantum Heisenberg model whose classical counterpart 
is given by the Hamiltonian in Eq. (1). In that case, one 
simply obtains 1/r = J/h. From Eq. ( |l2| ) we see that 
the total spin S is conserved during the dynamical de- 
velopment, and is thus constant both in magnitude and 
direction. Equations ( |l0| ) and (|ll]) are easily interpreted 
as describing the precession of the vectors S13 and S24 
about the constant vector S, keeping their lengths in- 
variant. We note that Eqs. (S) and (0) show that each 
individual spin executes a more complicated dynamics, 
precessing about the particular S13 or S24 that describes 
the sum of its near-neighbor spins, which is itself precess- 
ing about the constant S. From well-known examples of 
rigid-body dynamics, we thus expect that the motion of 
the individual spin vectors will feature two frequencies, 
one for precession about S, and the other for precession 
about either S ± S24 or S ± S13, respectively. 





x.s 4 


S/ S24; 













^24 + C*24 — $24- (16) 

We also have the relations 

^24 = ~Al3 (17) 

and 

C 13 + C 24 = S. (18) 
By combining Eqs. ([l5]) and (|l8|), we obtain 

o2 1 o2 c2 

C 13 = S + g 24 (19) 

and 

c-2 , o2 _ o2 

C 24 = + f s V (20) 

We now determine the individual spin vectors S^. Be- 
cause the four equations in Eqs. (||) and (g) have the 
same general structure, it suffices to focus on just one 
of them, say S2(t). We write S2 in terms of its compo- 
nents, 52s, S^x, and S2y, and make use of the standard 
Fourier transform S2i{t) = f 2 exp(zwt)5 2 i(w). We also 
let ^2± = i iS2y, and u>± = uj±S/t. We then obtain 

ojS 2s (u) = ^ [S 2 +(oJ-) - S 2 -(uj + )} (21) 
and 



FIG. 2. Sketch of the spin vectors. S, S13 and S24 all lie 
in a common plane where S13 and S24 each precess about S. 
The unit vectors Si and S3 initially point an arbitrary angle 
<?!>io out of the plane, and precess about S13. Similarly, the 
unit vectors S2 and S4 initially point an arbitrary angle (j>20 
out of the plane, and precess about S24. 

The conservation of S enables us to solve Eqs. ( |l(i| ) 
and ([ll]) exactly. Since S = S13 + S24 at all times, it is 
convenient to describe the motion in terms of the plane 
containing the three vectors S, S13, and S24, as pictured 
in Fig. 2. We thus have 

S 24 (£) = C 24 s + v4 2 4[xcos(st/r) - ysm(st/r)] (13) 
and 

Si 3 (i) = C13S + ^i 3 [xcos(st/r) - ysin(st/r)], (14) 

where s is a unit vector parallel to S, x and y are unit 
vectors normal to S satisfying x x y = s, thus completing 
the orthonormal basis set, and the constants Aij and Cy 
satisfy 



A-13 + C 13 



'13 



(15) 



and 



w5 2 ±(o;) = t— 52±H±— 5 2 ,(w±). (22) 

T T 

Solving for S2±{uj), and then replacing to by uj^, we have 



S2±(U T ) 



C24 ± LOT 



(23) 



where we have employed Eq. (^8[). Solving for £2s(w), 
we then find 



LO 



■(5 2 2 4 -(^r) 2 )52 S M -0. (24) 



This implies that 52s (w) has components from three fre- 
quencies only, lu = 0, ±<S24/t. Thus, we write 



S 2b (w) = 2wS2so5(lv) + 7rA5 2s0 ^ 20 5{uo - S 2 a/t) 



+e- l ^°5{io + S 24 /r) 



or in real time, 



(25) 



(26) 



where the constants <S2sO and A^so will be determined 
below, and </>20 is the arbitrary angle at which S2 initially 
makes with the Fig. 2. We note, however, that S 4 must 
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make the same initial angle 4>4o — 4> 2 q with this plane, 
since both spins have unit length, and their sum S24 is 
contained within that plane. Analogously, Si and S3 
both make the arbitrary initial angle <^>io with that plane. 
Combining Eqs. (|2^) and (|2^), we may solve for 



S 2 ±{uj) 



2ttA 2 4S 2s0 S(uj ± S/t) 



+e 



C24 

Tl4>20 S[u;±(S + S 2i )/ 



± lc f> 20 S[uj±(S - S 24 )/t] 

C24 — S 2 4 



C 



24 



'24 



(27) 



Inverting the Fourier transform, and using S 2x = (S2+ 
S2-)/2, S 2y = {S 2+ - S 2 -)/(2i), we have 



S 2x (t) 



24 



S 2s0 cos(St/r) 



'24 

A 2 4AS 2s q 

2[S 2 4 + C 2 4] 

A 2 4AS 2s o 
2[S 2 4 — C24] 



C0S[(5 + S 2A )t/T + 020 

COSTS' - S 2i )t/r 



20 



(28) 



and 



+ 



A 2 4S 2s q 

C24 
A 2 4AS 2s o 

2[S 2 4 + C24] 
A 2 4AS 2s o 
2[S 2 4 — C24] 



sm(St/r) 

sin[(5 + S 24 )t/T - 
sin[(5 - S M )t/T 



»20j 



»20j 



(29) 



Thus, as pictured in Fig. 2, S2, and correspondingly, S4, 
precesses about S24, which itself precesses about S. As 
in Fig. 2, Si and S3 each precesses about S13, which 
itself precesses about S. 

To evaluate the amplitudes S 2s q amd AS 2s q, we use 
the fact that S'l(t) = 1 to obtain 



<-,2 
°24 



2M) 



(AS- 



2s0) 



r 2 

°24 



A 2 
A 24 



= 1, 



(30) 



which is independent of 020- This equation provides a 
constraint upon the two amplitudes S2s0 and AS 2s q. We 
may then immediately write down the analogous equa- 
tion for the amplitudes S4 S q and AS4 S o appearing in the 
analogous expression for S4(i), 



5. 
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( ^ll 



(AS 4s0 ) 2 



A 2 4 



= 1. 



(31) 



Since S4(t) = S 2 4(t) — S 2 (t), where S 2 4(i) is given in Eq. 
(H), the expressions for S 2 (i) given by Eqs. ©, (Q, 
and (29), and the analogous ones for S4(t) (with S 2s q 
and AS 2s o replaced by S4 S o and AS^o, respectively) are 
consistent, provided that 



S4 



C 2 4 — S 2 



and 



ASisO = - AS; 



2s0- 



Substituting Eqs. @ and ([331) into Eq. ([ 
tracting the results from Eq. fl30|), we have 



(32) 
(33) 

ilh, and sub- 



S 2s — S4s0 — C24/2. 
Then, from Eq. (|3^), we find 

AS 2s0 = -AS is0 = 4^[1 - SlJA] 1 / 2 . 
024 



(34) 



(35) 



In Eq. (p5|), we have made the arbitrary choice of as- 
signing the positive sign to AS 2s q, but that does not 
affect any of the results. Thus, we have now completely 
determined the dynamics of S 2 (t) and S4(£), except for 
the arbitrary phase 4> 2 q representing the angle that S2(0) 
makes with the plane containing S, S24, and S13. Simi- 
larly, S x (i) is obtained from Eqs. ©, (||), (f§), @), 
and ( |35| ) by replacing A24, C24, S'24, and (f> 2 o with A13, 
Ci3j S13, and 4>iq, respectively. Ss(t) is then obtained 
from Si(i) in the same way as S4(t) was obtained from 

s 2 (t). 



III. TIME CORRELATION FUNCTIONS 

In this section, we utilize the exact results for the dy- 
namics of the four spin vectors derived in the previous 
section to obtain analytical formulas for the three dis- 
tinct time correlation functions. 



A. General Results 

There are three inequivalent correlation functions, 
which we denote by C 22 {t) = (S 2 {t) ■ S 2 (0)), C 12 (t) = 
(S 1 (t)-S 2 (0)), and C 24 (i) = (S 2 (t) • S 4 (0)>, where (...) = 
Tr[exp(— f3H) . . ]/Z. These are the spin-spin autocorre- 
lation function, the near-neighbor spin-spin correlation 
function, and the next-nearest-neighbor spin-spin corre- 
lation function, respectively. In evaluating these func- 
tions, we must average over the initial conditions, which 
means not only the averages over S, £24, and S13, but 
also over the initial angles 0io and cf> 2 Q , which are present 
in Eq. ([|) in the integrations over the solid angles fli and 
2 . We note that all of the correlation functions depend 
upon the temperature through the parameter a, but to 
keep the notation simple, we suppress that dependence. 
The simplest of these Cy(t) is C\ 2 (t) : for which the inde- 
pendent averages over 4>w and 4> 2 q greatly simplify the 
final expression. We find 



4 



Cia(t) = -{C 13 C 2 4 + A 13 A 24 cos(5i/r)) (36) 



-= dx dii 

*-(x 2 ~y 2 ) 2 



x+y 



/ 2 2 2\ 

sds e a(s - x - y ) 



x-y\ 



V 



As 2 

[s 4 + (x 2 -y 2 ) 2 ] 
4s 2 



) cos(si/r) , (37) 



where a = /3J/2, Z is given by Eq. (Q), and we have 
replaced 5 1 , 5i3, and 624 by s, x, and y, respectively. 
After some algebra, the autocorrelation function is found 
to be 



C 22 (t) =Io + h(t) + h(t) + h(t), 
where 

-fo = t(C| 4 ), 
7 1 (t) = i(A 2 4 cos(5t/r)), 

/ 2 W = 5(J i [l-5 2 2 4 /4]cos(5 2 4t/T)) ) 



(38) 

(39) 
(40) 
(41) 
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and 

hit) 



1, [1 -SIJ4] 
2^ <? 2 



([Cf 4 + Sf 4 ] x 



x cos(St/r) cos(S24t/r) 
+2C* 24 S*24 sin(5t/r) sin^t/r)) ). 

Using the same notation as in Eq. (|37|), we have 

1 Z" 2 Z" 2 /" a, +J/ 
Jo = -= / dx I dy sds x 

AZ Jo Jo J\x-y\ 



(42) 



xe 



«( s 2 -x 2 -j/ 2 ) (s 2 -^ 2 + y 2 ) s 



As 2 



1 



Ji(t) = -4= / dx I dy I sdse a( - s 

4Z Jo Jo J\x-y\ 
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x cos(si/r) (y 2 



(s 2 - x 2 + y 2 ) 2 



4s 2 



(43) 



(44) 



hit) 



7= dx dy 
£6 Jo Jo 



(1 - y 2 /4) 
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cos(yt/r) x 



r x+y 

x / sdse a( - s 
■ ! " 

„2 ™2 



(s 2 - x 2 + y 2 ) 2 
4? 



x y 



(45) 



and 

hit) 



1 

2^7 



2 /~ 2 (i-2/ 2 /4) r +y 

g?x / c?y t. / sds x 

y J\*-v\ 



n 



f 2 2 2\ 



(s 2 - x 2 + y 2 ) 2 
4s 2 



+ 



x cos(yt/r) cos(st/r) 
y(s 2 -x 2 + y 2 ) 



(st/r) sm(yt/r) 



Finally, C 2 4(t) may be found from 

<S 24 (i) -S 24 (0)) = 2C 22 (i) + 2C 24 (t), 



(46) 



(47) 



or 



Cuit) = -C22W + i(C 2 4 + A 2 24 cos(St/r)> 

= 2I + 2/ 1 (i)-C 22 (t) 

= I + I 1 (t)-I 2 (t)-I 3 (t). 



(48) 



By interchanging the integrations over x and y, it is 
easy to see that (Ci3C 2 4) = (S 2 /2 — Cf 4 ). Thus, we note 
that Ci 2 (i) is equivalent to 



CMt) = (S 2 )/8-I -h(t). 



(49) 



We therefore remark that the (t) satisfy the conserva- 
tion law, 

C22W + C 24 (t) + 2C 12 (t) = (S(t) • S(0)>/4 = (S 2 )/4, 

(50) 

a temperature-dependent quantity. Hence, in the infinite 
time limit, two of the three correlation functions Cij(i) 
approach the same limit, 



but 



lim t ^ 00 C 22 (Z) = lim t ^ tX) C 24 (t) = I Q , 



limt^daCt) = (S 2 )/8-/ , 



(51) 



(52) 



since the other terms vanish due to the infinite number 
of oscillations of the integrand within the interval of in- 
tegration. This is essentially a consequence of angular 
momentum conservation. 122] 



B. Reduction to Quadrature at Infinite Temperature 

1. Analytic results at infinite temperature and time 

In the limits t, T — > 00, we can evaluate the Cij(t) 
analytically. From Eqs. |58|), @, and @, we note 
that these three functions are all given by the T — > 00 
limit of Io and the Ii(t) for i = 1,2, and 3. We first 
consider the simplest of these, Jo, which gives the t — > 00 
limit. We find, 



lim C 22 (<) = lim C 24 (i) = \ + S 4 



(53) 



and 
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lira Cu{t) = - —Si, 



where 



Si = — In 2 
45 



11 
180 



0.062115. 



(54) 



(55) 



At first sight, one might have intuitively expected 
that the three Cij(t) should be equal to each other as 
t, T — > oo, and since lim T ^ 00 (S' 2 }/4 = 1, Eq. would 
require each of them to equal 1/4. This expectation is 
in fact the result predicted by conventional diffusive spin 
dynamics in the infinite temperature limit. p3f Moreover, 
for finite times that formalism predicts that all of the cor- 
relation functions depart from their common infinite-time 
limit by terms that decay exponentially to zero. However, 
our present rigorous results, Eqs. (53) - (pq) , as well as 
( |59| ) - ([H]) in the following, show that these expecta- 
tions are without foundation. Similar findings apply for 
the following simpler systems: the classical dimer, equi- 
lateral triangle, and regular tetrahedron, p3| , for which 
the exact time correlation functions are derived as one- 
dimensional integrals for all times and temperatures. 



dsh,3(s) cos(st* 



(58) 



where t* = t/r, and analytic forms for the /i(s), gi{s) 
and /i3(s) are listed in the Appendix. The infinite tem- 
perature correlation functions £22(^)7 C24(£), an< ^ ^12 ft) 
are then simply found using Eqs. (|38|), (fig), and (|4S|), 
respectively. Thus, we have reduced these infinite tem- 
perature correlation functions to quadrature. They are 
shown for < t/r < 10 in Fig. 3a. 



pj/2 = 



C12 

- c 22 (a) 
. C24 




t/x 



2. One- dimensional integral representations 1 

In this subsection, we give one-dimensional integr 
representations for the three time correlation functions 
infinite temperature. One important advantage of 
reduced forms is that they allow us to easily derive 
alytical formulas for the leading corrections to the Ion 
time asymptotic values for finite times of each of thej g 
relation functions. Another important advantages 
it becomes possible to obtain extremely accuratcoiumc r 
ical values for the infinite temperature correlation 
tions for all times. By comparison, for finite tempera&mie 
accurate numerical evaluation of the three-dimensional 
integrals in Eqs. (^3j) - ( ^6| ) becomes a major challengj 

We have found that the three functions Ii(t) may be 



written as 



0.2 



lim h(t) - 

1 — >oo 



lim I 2 (t) - 

1 — »oo 



and 



lim I 3 (t) 

1 — >oo 



ds fi{s) cos(si*) 



+ J ds gi(s) cos(st*), 
2 

dsf2{s) cos(si*), 



dsfs(s) cos(st*) 



+ J dsg 3 (s) cos(st*) 




(57) 



FIG. 3. (a) Plot of Ci 2 (t) (solid), C 24 (i) (dotted), and 
Cizit) (dashed) versus t/r in the infinite temperature limit 
a — 0. (b) Comparison of the numerical results of Muller 
[17] (dashed) with our exact results (solid) for the infinite 
temperature (a = 0) autocorrelation function C22(t). 

We remark that the infinite temperature autocorre- 
lation function C22{t) was obtained previously using a 
purely numerical procedure. |17| In Fig. 3b, we have 
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compared those published results with our exact formula 
at infinite temperature. Although there was some dis- 
tortion in the axes in the published figure, using a pure 
rotation to account for this distortion led to the excellent 
agreement between the numerical and exact results. 

From Fig. 3a, the autocorrelation function C22W de- 
creases from its initial value £22(0) = 1, then undershoots 
its aymptotic limit i + 84, and approaches this limit by 
oscillating about it for a rather long time. On the other 
hand, spins of different sites are initially uncorrelated at 
infinite temperature, Cn{Q) = €24(0) = 0. At later limes, 
these functions both overshoot their respective asymp- 
totic limits \ — 84 and \ + 64, and then oscillate about 
them. The oscillations of Cx2(t) decay so rapidly that 
they are barely discernible in this figure. On the other 
hand the oscillations of C24W are of the same amplitude 
and persist as long as do those of C2i{t), and are likewise 
easily seen in this figure. In addition, after C22W and 
C24(t) first become equal to each other, they braid about 
each other in their approaches to the same aymptotic 
limit. 

In order to see more clearly how this occurs, we have 
found analytic expressions for the leading behaviors of 
the correlation functions for long times, t » r. We 
first consider Cyi{t). In this case, besides the constant 
Io, we only need to evaluate Ii(t). To do so, we inte- 
grate both terms in Eq. ( |56| ) by parts, treating cos(st*) 
as the variable to be integrated, and /i(s) and g\ (s) as 
the variables to be differentiated. From the results in the 
Appendix, it is seen that /i(s) and gi{s) as well as their 
first three derivatives are continuous at s = 2. In addi- 
tion, since the relevant integration endpoint values and 
derivatives at s = and s = 4 also make no contribu- 
tion through third order in the repeated integrations by 
parts, the leading contribution to the aymptotic behavior 
arise from the non-vanishing f 1 (0) and g 1 (4). The final 
result for the leading behavior is given by 



Jim Ci 2 (t) 



*4+4^4[|-cos(4f)]. CO, 



On the other hand, C22W and Cviif) at infinite tem- 
perature also depend upon I2 (t) and / 3 (t) . Again, we in- 
tegrate by parts in a similar fashion, treating /2(s), /a(s), 
53 (s), and /i 3 (s) as the variables to be differentiated. The 
leading non- vanishing contributions to C22W and C24W 
from these integrations by parts are both of second or- 
der. For l2(t), the leading non- vanishing contribution 
comes from the non- vanishing f 2 (0) and j ' 2 (2), the latter 
of which is a non-trivial number. For the leading 

non- vanishing contribution arises from / 3 (0), / 3 (2) and 
<7 3 (2). Although both / 3 (s) and g 3 (s) have non-trivial 
values and derivatives at their matching point ,s = 2, 
the difference between their derivatives is a trivial, but 
non- vanishing value. In addition, the functions f^s) and 
<? 3 (s) both have non-trivial values and derivatives at their 



matching point s — 4, but these values and derivatives 
are equal, and thus their contribution in second order to 
the integration by parts vanishes. We thus obtain the 
long-time behaviors at infinite temperature, 

Km C 22 (t)^ ^ +64 



T — >co 
t>> 



and 



lim C24W 



T — > 

t>> 



+ 60f^ [5_(29 + 81n2)c ° s(2r)] ' (60) 



1 . 
- + 84 



1 



r [5-(29 + 81ri2)cos(2t*)]. (61) 



60t* 

We note that C\2{t) decays much more rapidly (oc 
1/i* 4 ) to its constant long-time limit than do either C22W 
orC24(£) (oc 1/i* 2 ). The long-time braiding of these func- 
tions about each other arises from the opposite signs 
of their oscillatory terms. Furthermore, at long times, 
C12 it) oscillates with twice the frequency of the long-time 
oscillations of C22W and C24(t). 



C. Results for finite temperatures 

At finite T, we evaluate Io{a) and (S 2 )(a) numerically. 
For I a I < 1, breaking each integral into 100 intervals is 
sufficient to obtain 0.1% accuracy. Note that this means 
there are 10 6 integration intervals overall. However, for 
low T (| a I > 1), the number of intervals necessary to ob- 
tain that degree of accuracy increases. At \a\ — 10, one 
needs to break up each integration domain into 400 inter- 
vals, for instance. In Fig. 4, we have plotted the infinite- 
time limit of the spin-spin correlation functions Iq (a) and 
(<S 2 )/8 — Io(a) for both the ferromagnetic (FM) and an- 
tiferromagnetic (AFM) cases. As a — > 0, one obtains the 
analytic limits given by Eqs. (|5^) and (|54|). However, in 
the low T limit \a\ -> 00, both I and (stf/S - 1 -> 0(1) 
for the AFM (FM) case, respectively. This just tells us 
that at T = 0, all of the spins are aligned in the FM case, 
and in the AFM case, their sum is 0. We note that for 
a j << 1, Io(a) obeys the inversion symmetry, equiva- 
lent to dlo/da exists. We note that for the AFM 

case, lim^oo Ci2(<) is negative for a < —0.71. This just 
reflects the fact that for the antiferromagnetic ring, the 
spins on neighboring sites are anticorrelated at long times 
and for temperatures that are not too large. 
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|a| 

FIG. 4. Plot of 7 (a) = lim t _ 00 C 2 2 (*) = lim t _ >00 (724(*) 
and {S 2 }/8 — To (a) = lim t _ 00 Ci2(i), as a function of \a\, for 
the FM (a > 0) and AFM (a < 0) cases. 

At finite T we also may evaluate the Ii (t) numerically 
from the triple integral forms, Eqs. ( (44| ) - (fl6|). As for Iq, 
we break each of the three integrals into N intervals. At 
the lowest T values considered (a = —20), it is necessary 
to take N > 1000 to achieve sufficient accuracy. In Figs. 
5-7, we have plotted the Cij(t) for the FM case with a — 
0.5,2, and 10, and compared with the analytic results 
for a = 0. In each of these figures, C\2{t) decays to 
the equilibrium value (S 2 )/8 — Io(a) more rapidly than 
and C2i(t) decay to their mutual equilibrium value 
Io{a), while oscillating for a few periods about the latter. 
As T decreases, all of the Cy(0) increase monotonically, 
approaching unity as T — > 0. In addition, the oscillations 
persist to much longer times. Also, as seen in Figs. 5 
and 8, as T decreases, C\2(t) oscillates for an increasing 
amount of time, and C22W and C24.it) oscillate about it 
for an even longer period of time. 
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FIG. 5. Plots of & 2 (t) versus t/r for the a = /3J/2 FM 
cases 0, 0.5, 2, and 10, as indicated. 
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FIG. 6. Plots of C 22 {t) versus t/r for the a = (3 J/2 FM 
cases 0, 0.5, 2, and 10, as indicated. 
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t/T 

FIG. 7. Plots of C 2 4.(i) versus t/r for the a = /3J/2 FM 
cases 0, 0.5, 2, and 10, as indicated. 

At lower T, the amplitudes of the oscillations eventu- 
ally reach a maximum, so that the oscillations in C\2(t) 
for a — 0.5, 2, 10 are distinctly noticeable. As T — > 0, 
the lifetimes of the oscillations appear to diverge, but 
their amplitudes become vanishingly small. At a = 10, 
we have shown the behaviors of C\2{t) (solid), C22W 
(dashed), and C^iif) (dotted) together in Fig. 8 for the 
extended time domain < t/r < 30. Throughout this 
domain, the decay of the oscillations in all three correla- 
tion functions is small but discernible. However, careful 
inspection of the oscillating waveforms reveals that C12 (t) 
oscillates with twice the frequency of the other two, con- 
tinuing the pattern that we have already seen for infinite 
temperature. Note that Cxzif) appears to oscillate nearly 
as a simple cosine function, but Cmfo) and C2i{t) have a 
more complicated oscillatory behavior, with a fundamen- 
tal frequency that is one-half that of C12 (t) , and they are 
almost completely out of phase with respect to one an- 
other. 
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FIG. 8. Plot of & 2 (t) (solid), C 24 (i) (dotted), and C 22 (t) FIG. 10. Plot of C 12 (t) (solid), C 24 (t) (dotted), and C 22 (t) 

(dashed) versus t/r for the low temperature FM case (dashed) versus t/r for the AFM case a = —2. 
a = /3J/2 = 10. Note that C 22 (0) = 1 and C 24 (0) » 0.95. 



The corresponding results for the AMF case are shown 
in Figs. 9-11, for which a = -0.5, -2, and - 20, respec- 
tively. The last of these, a = —20, took weeks of compu- 
tational time to obtain sufficient accuracy. In these cases, 
we presented the Ci2(t), C22(i), and C2i(t) data as solid, 
dashed, and dotted curves, respectively. We note that 
as the temperature is lowered, £12(0) decreases towards 
the value -1, which would correspond to perfect AFM be- 
havior. However, C\2(t) then increases with t, reaches a 
maximum, and then decreases to the asymptotic , infinite 
time limit. In addition, as T is lowered, C24(0) increases 
towards +1, approaching £22(0). Then, at some time t\, 
£22(^1) nrs t equals C2i(ti), an d thereafter, the two func- 
tions are braided about each other. The braiding oscilla- 
tions decrease in amplitude as T is decreased, so that the 
overall dj (t) all approach non-oscillatory uniform curves 
as T -> 0. 
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FIG. 11. Plot of Ci 2 (t) (solid), C 24 (t) (dotted), and C 22 (t) 
(dashed) versus t/r for the very low temperature AFM case 
a = -20. 
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FIG. 9. Plot of & 2 (t) (solid), C 24 (t) (dotted), and C 22 (t) 
(dashed) versus t/r for the AFM case a = —0.5. 



From numerical simulation studies of more compli- 
cated mesoscopic classical systems, it has been suggested 
that the low temperature AFM autocorrelation function 
should scale, approaching uniform functions of tT 1 / 2 . ]24j] 
To investigate whether such a scenario holds for this ex- 
actly solved four-spin system, we first plotted the AFM 
autocorrelation function at the low T values we consid- 
ered. This is shown in Fig. 12. Although there are 
oscillations that persist to increasing times as T is low- 
ered, the overall shape of the curves does not change its 
shape qualitatively, suggesting that C^if) might indeed 
scale as a single function of tT 1 / 2 as T — > 0. In fact, 
this behavior has been established by analytical means 
for the simpler cases of the classical dimer, equilateral 
triangle, and regular tetrahedron. [2l|j2^ | 
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FIG. 12. Plot of C 22 (t) for the a = /3J/2 AFM cases -0.5, 
-1, -2, -5, -10, and -20, correspondingly from top to bottom 
at large t/r, as indicated. 

In Fig. 13, wc therefore plotted C22W versus 
t/[T|/3J/2| 1 / 2 ], which is proportional to tT 1 / 2 , to check 
this notion quantitatively. Indeed, the curves do scale, 
except for the braiding oscillations, which are decreas- 
ing in magnitude as T decreases. Thus, curves for the 
lowest two temperatures, a = —10 and a = —20, would 
nearly fall on top of each other if the oscillations were 
not present. Similar low temperature scaling behavior of 
C2i(t) is shown in Fig. 14, which also includes the braid- 
ing oscillations. C\2{t) exhibits a clearer example of the 
scaling, as shown in Fig. 15, since it does not contain any 
braiding oscillations. The major deviation from scaling 
occurs at very short times, although the differences be- 
tween the curves at a = —10 and a = —20 are not so 
large there. Since £12(0) = — 1 in the zero temperature 
limit, this deviation from scaling probably arises from the 
fact that the T = limit has not yet been reached for 
such short times. 
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FIG. 13. Plot of C2-2.it) as a function of the scaled time 
(t/r)|/3J/2| 1/2 , for the a = (3J/2 AFM cases -0.5, -1, -2, -5, 
-10, and -20, correspondingly from top to bottom at large t/r, 
as indicated. 
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FIG. 14. Plot of C2i{i) as a function of the scaled time 
(t/r)|/3J/2| 1/2 , for the a = (3J/2 AFM cases -0.5, -1, -2, -5, 
-10, and -20, correspondingly from top to bottom at large t/r, 
as indicated. 
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FIG. 15. Plot of Ci2(i) as a function of the scaled time 
(t/r)|/3J/2| 1/2 , for the a = (3J/2 AFM cases -0.5, -1, -2, -5, 
-10, and -20, correspondingly from top to bottom at small 
t/r, as indicated. 

At higher temperatures, the scaling property gradually 
breaks down. This is clearly seen for a = —0.5 in Figs. 
14 at 15, for which the short-time values of C2i{t) and 
C\2{t) deviate greatly from the values (1 and -1) obtained 
respectively in the low temperature limit a — > — 00. For 
C\2{t), the deviations are also rather large at long times. 
However, for all three correlation functions, even at a = 
—0.5, the positions of the dip for €22^) and C2i{t) and 
the peak for C\2(t) still scale. 

D. Fourier Transforms 

We now evaluate the Fourier transforms Ii(uj) of the 
Ii(t). Since C22W an d C2i(t) both approch the constant 
Io i= as t — > 00, and C\2{t) approaches {S 2 )/8 — Io i= 
in the same limit, the Iq ^ or (S 2 )/8 — Iq ^ present in 
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the respective Cij (t) give rise to delta functions in terms 
of the angular frequency oj, in Cij (00), equal to either 
2ttI 5(oj) or 2n((s 2 )/8 — Iq)S(oj), which can be written 
down by inspection from Eq. (Jl3|). We shall therefore 
evaluate the Fourier transform of the deviations 5Cij (t) = 
dj(t) - limt^ooCy ^), 



= 2ttI 6{lj) + Sdjiuj) 
for (i,j)=(2,2) and (2,4), and 

C 12 (lo) = 2vr((5 2 )/8 - I )6(oo) + <yCi 2 (w), 
where 



SCij(u;) 



dte- iujt 6Cij{t). 



(62) 



(63) 



(64) 



Causality requires that 5Cij(—oo) = SCij(uj). Thus, it 
suffices to evaluate the 5Cij(co) for 00 > 0. For simplicity, 
we first evaluate 12(00), and let a) = wt. Since the only 
time-dependence within the expression for I2 (t) appears 
in the factor cos(yt/r), Fourier transformation replaces 
this factor with [8(00 — y/r) + S(oo + y/r)]/2, which can 
be written as (r/2)[6(y — Co) + S(y + Co)]. The second 5- 
function does not contribute to the Fourier transform for 
00 > 0, so we obtain for oj > 0, 

U.) = ^e(2-o»(l-a)V4) [> ^ f ^ z x 



8C0Z 



J -1 

2 ,~,2\2 



{ 9 -9 „ - (X — 0) ) 

xlx +ui — 2xzuj - 



x 2 + Co 2 + 2xzuj 



(65) 



where Q(x) is the Heaviside step function, and we have 
made the change of variables s 2 — * x 2 + Co 2 + 2xzCo for 
ease of computation. 

We now evaluate Ix(oo). As in the expression for 12(00), 
Fourier transformation of the factor cos(si/r) replaces it 
with (t/2)[5(s — u>) + S(s + &)], and the second term 
does not contribute to the integrals for ui > 0. However, 
since < s < 4, there are now two regions of integration 
over the variables x and y. For 2 < ij < 4, the only 
region of integration is the interior of the isosceles triangle 
with sides obeying x — 2, y — 2, and x + y = Co, and 
corners at their intersections. For < Co < 2, the region 
of integration is the interior of the pentagon with sides 
obeying y — 2, x — 2, y — x = —Co, y + x = Co, and 
y = x + Co. This interior region is symmetric about the 
line y — x, and can be broken up into two regions of 
integration. The first region is the interior of a rectangle 
rotated 45° about the axis normal to the xy plane, with 
sides obeying y = x ± Co and x + y = 2 ± (2 — Co). The 
second region is the interior of the isosceles triangle with 
sides obeying y = 2, x = 2, and y = —x + 4 — Co. 



In the triangular integration regions, we maintain the 
integration variables x and y, keeping account of the in- 
tegration limits. However, in the rectangular integration 
region, it is convenient to perform a rotation of the axes 
by 45°, letting r = x — y, s = x + y, and incorporating 
the Jacobian, which replaces the differential integration 
area dxdy with drds/2. We thus have 

7 1 M = 7rrQ( "- 2) 5 4 -^ f dxf dyx 

I&LOZ Jco-2 Ju-x 

xe a ^ 2 - x2 - y2 1{-Co 4 

+2Co 2 (x 2 +y 2 )-(x 2 - y 2 ) 2 } 

^ 2 -^ f 2 dxT dyx 

2—ui J 4 — x—u 



I6C0Z 



-Co 4 



+2Co 2 (x 2 +y 2 )-(x 2 -y 2 ) 2 } 

+ " T6(2 1" ) [ Q dr r^dse^ 2 -^ 2 ^ x 
I6C0Z Jo Jcj 

x[-Cb 4 + Cb 2 {r 2 + s 2 )-r 2 s 2 }. (66) 

Last, but by no means least, wc evaluate 13(00). This is 
easiest to do if we first use the trigonometric relations to 
rewrite Eq. ( |46| ) in terms of cos[(y + s)t/r] and cos[(y — 
s)t/r]. Then, after Fourier transformation, we obtain, 



h(oS) = — = / dx dy 
IbZ Jo Jo 



(l-y 2 /4) f x+v ds 



xe" 1 " - * } (j(y + s) 2 - x 2 ] 2 [5(oj + y - s) 

+6(Co-y + s)} 
+ [(y - s) 2 - x 2 } 2 {5(oj + y + s) 

+S(oj — y — s)] 



(67) 



For oj > 0, the term containing 8(uj+y+s) vanishes, so we 
are left with three terms, which we examine separately. 
We denote them I^n, where n = 1,2,3, corresponding 
to the order in which the remaining <5-functions appear 
in Eq. (|67|). The first integral, /31, is subject to the 
constraints \x — y\ <Co + y<x + y. The second (Co + y < 
x + y) of these two constraints implies Co < x, which also 
implies < Co < 2. If y < x, the first constraint implies 
y > (x — Co)/2. On the other hand, if x > y, there is 
no additional constraint on the integration region, other 
than x > Co. Thus, the combined integration region of 
this integral is the interior of an irregular quadrangle with 
sides obeying x = Co, y = 2, x = 2, and y — (x — Co) /2. 
Hence, we may now write l3i(to) by inspection, 



hi(oo) 



7tt6(2 - Co) 



dx 



16Z jcj J{ x 

3 a(Cj 2 +2uy-x 2 )^2y 



s)/2 y(u + y) 



Co) 2 



„2l2 



(68) 



The integration regime of the second integral, 132(00) 
is subject to the constraints \x — y\ < y — Co < x + y. 
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The second "constraint" makes no restrictions for Co > 0. 
However, the first constraint can be rewritten as Co < 
y — \x — y\, and since x and y are both between and 
2, this then implies < Co < 2. For x > y, the first 
constraint further requires y > (Co + x)/2. On the other 
hand, for y > x, the first constraint requires x > Co. 
Thus, the integration region is the interior of the irregular 
quadrangle with sides obeying x = Co, y = 2, x = 2, and 
y = (Co + x)/2. We therefore write 



h2 (oj) 



7rr6(2 



Z Jcj J(Co+ 



(1 - j/ 2 /4) 



xe 



16Z 

a(u) 2 — 2yuj — x 2 ) 



d y- 2i ~\ 

x)/2 y 2 (y-u) 

~.\2 2l2 



[[2y-CbY-x 2 ] 



(69) 



Finally, the integration region of i33(w) is subject to 
the constraints \x — y\ < Co — y < x + y. The second 
constraint implies that Co can be as large as 6, and also 
that y > (Co — x)/2. For x > y, the first constraint implies 
x < u>, whereas for y > x, it implies y < (Co + x)/2. 
For < Co < 2, these constraints restrict the integration 
region to the interior of the triangle with sides obeying 
y = (Co — x)/2, y = (Co + x)/2, and x = Co. For 4 < Co < 6, 
the integration region is the interior of the triangle with 
sides obeying y = 2, x = 2, and y = (Co — x)/2. In the 
intermediate regime 2 < Co < 4, the integration region is 
the interior of the irregular quadrangle with sides obeying 
y = (Co — x)/2, y = (Co + x)/2, y = 2, and x = 2. We 
break this integration region up into two parts. One of 
these parts is the interior of the isosceles triangle with 
sides obeying y = (Co — x)/2, y = (Co + x)/2, and x — 
4 — Co. The second region is the interior of the irregular 
quadrangle with sides obeying x = 4 — Co, y = 2, x = 2, 
and y = (Co — x)/2. Altogether, we write hz(oo) as 



hz(w) 



nTe(2-Co) [* r(-+-)/ 2 

ax I dyf(x,y,u) 



16Z 



(w-x)/2 



htQ(Co - 2)0(4 - Co) 
+ 161 X 

4-tD MQ+x)/2 

dx / dy 

i! J(Q-x)/2 



f dx f 



dy 



4-lD J(di-x)/2 

tttO(Co - 4)0(6 - Co) 

16l X 

r-2 r 2 



f{x,y,Cj) 



f dx f dyf(x,y,u>), (70) 

JQ-A J(Q-x)/2 



'(O-x)/ 

where 

f(x,y,Co) = ^fI^e^-^-^ 
y 2 (uj - y) 

x[(2y-Co) 2 -x 2 } 2 . 



(71) 



In Figs. 16 - 20, we present our results for the SCij(u)), 
plotted as functions of lot. In the figure labels, we drop 



the tildes for clarity. Note that these functions are rigor- 
ously zero for lot > 6, but they are so small for lot > 4 
that they are indistinguishable from zero for lot > 4.1 
in plots with < lot < 6. In each case, 5Ci2(lo) are the 
solid curves, 8C2a(lo) are the dotted curves, and 8C22(u) 
are the dashed curves. At infinite temperature, these are 
obtained from our exact formulae for the fi(s), gi(s) and 
/i3(s) listed in the Appendix by simply letting s — > lot, 
and multiplying the overall results by ir. For example, 
5C\2(lo) = — 7t[/i(o;t) + gi(ioT)]. The exact expressions 
for the three 8Cij(oo) are plotted in Fig. 16. As a check 
on our computations, we also obtained these results by 
numerically evaluating the SCij(t) from the double inte- 
grals in the infinite temperature limit, and the results 
were found to agree to within three significant figures. 
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FIG. 16. Plots of the exact Fourier transforms 8C\i(lo) 
(solid), 8C22(lj) (dashed), and 5C2a(lo) (dotted) of 
Cij(t) — limt^ooCV,^), as functions of lot, in the infinite tem- 
perature limit a — 0. 
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FIG. 17. Plots of the Fourier transforms 8C\2(lo) 
(solid), 8C22(lo) (dashed), and 8C2a(lo) (dotted) of 
Cij(t) — \im t ^ooCij(t), as functions of lot, for the FM case 
a = 0.5 
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FIG. 18. Plots of the Fourier transforms 8C\2{lo) 
(solid), <5C22 (i^) (dashed), and SC 2 4 {lo) (dotted) of 
Cij(t) — \imt^ooC'ij(t), as functions of lot, for the FM cases 
a — 1,2,5, and 10. The peaks developing successively near 
to lot = 2, 4 with decreasing T are magnons. 
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FIG. 19. Plots of the Fourier transforms SC\2{lo) 
(solid), SC 2 2{^>) (dashed), and SC 2 4(lo) (dotted) of 
Cij(t) — lim t ^ooCij(t), as functions of lot, for the AFM cases 
a = —0.5 and -1. The bumps developing for a = — 1 at 
lot w 2 are magnons. 



In this and in subsequent figures, we also checked the 
accuracy of our analytic formulae by the zero-time sum 
rule, the J °° SCij (co)d,LO / V = 5Cij(0). It is seen that the 
functions all are continuous, and approach zero at large 
lo. SCi2(lo) < 0, and 5C22{oo) > 0, but 5C2a{lo) has re- 
gions of both signs. In Figs. 17 and 18, we present the 
data for the FM case, with a — 0.5 in Fig. 17, and 
a = 1,2,5, and 10 in Fig. 18. As the temperature is 
lowered, the peak in (5C 2 2(w) moves to higher lo values. 
Moreover, <5C22(^) and 5C2a{ U} ) approach each other in 
the region 2 < lot < 4, but approach the opposite of 
each other in the regime < lot < 2. Thus, in the regime 
2 < lot < 4 of Fig. 18, these dashed and dotted curves 
combine to give a curve that appears to be dash-dotted. 
In addition, as the temperature is lowered, the curves all 
develop into sharp peaks at lot w 2 and 4, which are 
asymmetric, dropping rapidly to zero at lot = 2 and 4, 
but having longer tails at lower lot values. More precisely, 
5Ci2(co) becomes a single peak at lot « 4, whereas the 
other two, <5C 2 2(u;) and 5C24(to) have identical peaks at 
lot w 4, which are opposite to that of 5C\2{lo). But, they 
also have larger and sharper peaks at lot « 2, which are 
opposite in sign to each other. These peaks at lot k, 2, 4 
arise from magnons. Hence, Fig. 18 provides a simple 
explanation of the detail of Cij(t) for a = 10 shown in 
Fig. 8. 



8 " f\ -10 

6 - ; \ 5Ci2 

/ A 

4-!/ \\ ----5C22 




0.5 1 1.5 2 2.5 



COT 

FIG. 20. Plots of the Fourier transforms 8C12 {co) 
(solid), 8C22(lo) (dashed), and 8C24(lo) (dotted) of 
Cij(t) — linit^oo CVj(i), as functions of lot, for the AFM cases 
a = —2, —5, and -10. The peaks for lot w 2 that successively 
sharpen with decreasing T are magnons. 

Curves for the AFM case are shown in Fig. 19 and 
20. In Fig. 19, we display the results for a = —0.5 and 
-1 together, and in Fig. 20, the results for a = —2, —5, 
and -10 are shown. As the temperature is lowered, the 
peak in ££22 (w) moves to lower frequency, resulting in 
the slowing down seen in the real time curves. It devel- 
ops into two peaks, a large one at low lot, and a small one 
at lot w 2, which is an antifcrromagnetic magnon. Note 
that the magnon width sharpens as T is lowered. Sur- 
prisingly, 5C24:{<-o) changes dramatically, mostly changing 
sign near lot « 2, but the small lo behavior increases to 
join the small lo behavior of 5C22{oo), and the opposite 
of it in the region lot w 2. These peaks at lot k, 2 arise 
from antiferromagnetic magnons. In addition, 8C\2(io) 
stays negative, and develops into a negative peak at low 
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lo which is opposite to that of 8C22 ( w ) and 8C2a(u))- 

Finally, in order to elucidate the nature of the slowing 
down as T is lowered, in Fig. 21 we plotted our AFM 
results for a = —0.5, —1, —2, —5, and -10 together, as a 
function of the scaled frequency ujt\PJ/2\ 1 ' 2 . We note 
that the position of the low frequency peak does indeed 
scale, but since the magnon appears at a fixed frequency, 
it does not scale. The magnon is distinctly visible for 
a = —2,-5, and —10 in this figure, as noted by the 
arrows. 
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FIG. 21. Plots of the Fourier transforms <5Ci2(w) 
(solid), <5(722 (w) (dashed), and 5624(0;) (dotted) of 
Cij{t) — \im t -*ooCij(t), as functions of the scaled frequency 
uit\PJ/2\ 1/2 , for the AFM cases a = -0.5,-1,-2,-5, and 
-10. The magnon peaks for cut « 2 do not scale, and are 
marked with arrows. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have presented the exact solution 
for the thermal equilibrium dynamics of four classical 
Heisenberg spins on a ring. For this system, there are 
three relevant time correlation functions, which are the 
auto- (C22), near-neighbor (C12), and the next-nearest- 
neighbor (C24) correlation functions, respectively. Using 
our results, we wrote exact expressions for these three 
functions in terms of triple integrals. At infinite tem- 
perature, we reduced these triple integrals to single inte- 
grals. We also obtained analytic expressions for the long- 
time, infinite temperature behavior of the three correla- 
tion functions. We found that the near-neighbor correla- 
tion function C\2(t) is strikingly different from the auto- 
and next-nearest-neighbor correlation functions C 2 2{t) 
and C 2 i{t). Not only does it approach a different finite 
value asymptotically, but it also oscillates with twice the 
frequency, and the oscillations decay in amplitude much 
more rapidly. Although the long-time asymptotic val- 
ues of these functions vary with temperature, the near- 
neighbor correlation function similarly differs from the 
other two functions for all temperatures. 



In addition, we were able to obtain the Fourier trans- 
forms of the deviations of the correlation functions from 
their infinite time aymptotic limits in terms of double in- 
tegrals. At infinite temperature, exact analytic forms for 
these Fourier transforms were obtained. As the tempera- 
ture is lowered, peaks in the Fourier transforms appear at 
lot = 2 for antifcrromagnetic coupling, and at lot = 2, 4 
for ferromagnetic coupling. These peaks sharpen up as 
the temperature is lowered. Although the origin of these 
peaks is purely classical, they correspond precisely to 
magnons, which are usually thought of as quantum me- 
chanical in origin. Here the magnons arise from standing 
waves, such as those on a violin string. These standing 
waves arise from the periodicity of the ring, as a com- 
bination of traveling waves moving both clockwise and 
counterclockwise in direction. For the antifcrromagnetic 
case, neighboring spins are opposite in direction at low 
temperatures, so that two full wavelengths fit into the 
ring. For the ferromagnetic case, one can have either two 
or four full wavelengths in the ring, the latter correspond- 
ing to every spin pointing in the same direction. 

In the quantum mechanical analogue, the magnon en- 
ergies are E n = 2J[1 — cos(fc„a)], where k n = 2m:/ L and 
L = 4a. We thus get E Q = 0, E x = 2J, and E 2 = 47. 
These latter two values correspond precisely to lot = 2, 4. 

In addition, for antiferromagnetic coupling, a second, 
much larger peak in the Fourier transform functions ap- 
pears at increasing lower frequency as the temperature is 
lowered. This second peak is found to scale with w/T 1 / 2 , 
corresponding to the iT 1 / 2 scaling of the low temperature 
antifcrromagnetic time correlation functions. 

It is interesting to compare the long-time asymptotic 
results of the correlation functions at infinite temperature 
with the known results for other Heisenberg rings. For 
a dimer, with two Heisenberg spins interacting via Eq. 
(1) with S3 = Si, the infinite temperature limit of the 
auto- and near-neighbor correlation functions were given 
for all t, [EBI 



Inn C llW 4- 1 + 2 o C °f^ 

T^oo w 2 2t* 2 



3sin(3f) 3[l-cos(2f)] 



2t* 3 

1 - lim Ci 2 (t), 

T — >oo 



-it* 



where we have used the conservation law, 
C u (t)+C 12 (t) = (S 2 )/2, 



(72) 
(73) 
(74) 

(75) 



and limy^ 0O (S' 2 )/2 = 1 for the dimer. 

For the case of three spins on a ring, it is not so triv- 
ial, but it is still much easier to evaluate the correlation 
functions at large times and infinite temperature than 
for the four-spin ring. The infinite temperature aymp- 
totic limit result for the autocorrelation function Cu(t) 
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was quoted previously, pq | and its complete derivation 
for all temperatures was given. p3| In this case, reduction 
to quadrature is rather simple, and one obtains a result 
for the three-spin ring 8Cn(t) analogous to Eq. ([56]). At 
infinite temperature, |23j 



lim SCn(t) 



dsf4(s) cos(st*) 



+ / dsgi(s) cos(st*), 



(76) 



where fi(s), (s), and their values and relevant deriva- 
tives at the integration endpoints are given in the Ap- 
pendix. Thus, for the three-spin ring, the leading terms 
for long times were shown to be, p3[ 



1 [sin(r)+sin(3t*)] 
km Cn{t) = - + 8 3 ^ . 

where 



(77) 



S 3 = ^-hi3- — w 0.147188. (78) 
40 10 ' 

In analogy with Eq. (|50|) , it is then easy to see that 



2C 12 (t)+C n (t) = (S 2 )/3, 
and hence that 



(79) 



since limT_ >00 (5' 2 )/3 = 1 for the three-spin ring. 

For N > 3 classical Heisenberg spins on a ring, it is 
straightforward to obtain the conservation equation for 
the correlation functions, 

(N + l)/2 

C n (t) + 2 ]T C ln (t) = {S 2 )/N, (81) 



71=2 



and 



N/2 



C n (<) + 2 ^ C iB (t) + C UN+2)/2 (t) = {S z )/N, (82) 

n=2 

for N odd and even, respectively. 

We note that for 1 < N < 4 at infinite temperature, 
{S 2 )/N = 1. More generally, at infinite temperature, 
the evaluation of (5 12 ) for an A^-spin ring maps onto that 
of the mean square displacement of a chain of length N 
during a random walk in three dimensions, and hence rig- 
orously lim T ^oc{S 2 )/N = 1 VN > 1. We then note that 
the long-time behavior of Cn(t) at infinite temperature 
arises primarily from the discontinuities in the derivatives 
of the functions fi(s), <?i(s), hi(s), etc. at the endpoints 
of the integration intervals. For N — 2, 4, the functions 
are finite and continuous over the integration regions, but 



at least one of their first derivatives is discontinuous at 
one or more of the integration endpoints. In addition, 
for the particular odd-spin ring with N = 3, both the 
functions and their first derivatives are continous at the 
integration endpoints, but the second derivatives are dis- 
continuous at the endpoints. The mathematical forms of 
the functions fi(s), ffi(s), hi(s), etc. and their various 
derivatives at the integration endpoints become increas- 
ingly complicated with increasing N for N < 4. 

It would be interesting to find out whether this "pat- 
tern" of matching the functions and their derivatives at 
the integration endpoints might be maintained for much 
larger N values. Thus, at least two possible scenarios 
that might develop from attempting to generalize our 
results to much larger N values arise. Regardless of 
whether N is even or odd, Cn(t) might behave for t* ~S> 1 
as 



Cn(t) 



Jj+6 N + J2n=Q a n COs(nt*)/r 2 , 



(83) 



which would occur if the functions were finite and contin- 
uous but with one or more discontinuous first derivatives, 
as for N = 2,4, or if they were finite but discontinous at 
one or more of the integration endpoints, respectively. Of 
course, behavior such as for N = 3 could also be obtained 
for higher N values, as well, as well as more complicated 
scenarios. We expect that limjv-«x> $n — 0, and that 
5 n decreases to zero faster than 1/N. We note that for 
N = 2, we have 82 = 0. 

Numerical simulation data for the N — 4,6,8,10,50 
autocorrelation functions were presented for < t* < 10 
by Miiller. [0 Those results show that the oscillations 
in Cn(t) for N > 4 are easily discernible out to t* = 
10. Hence, it appears unlikely that the autocorrelation 
functions would approach their asymptotic limits more 
rapidly than l/t* 2 . The high accuracy of the numerical 
results is exemplified by the excellent agreement between 
the result for N = 4 and our own exact analytical results, 
fL7[ as shown in Fig. 3b. 

The main difference between rings with even and odd 
numbers N of spins lies in the long-time behaviors of the 
two-spin time correlation functions. For both even 
and odd, the conservation law requires at least one of 
the Ci n for n ^ 1 to compensate for the leading long- 
time behavior of Cn(t). For odd N, we anticipate that 
all of the correlation functions will fall off with oscilla- 
tory corrections that have an amplitude of order (t*)~ m 
at large times, where m is likely to be a small natural 
number. We then raise the question for even N, as to 
which, if any, of the correlation functions approaches its 
asymptotic limit more rapidly than the autocorrelation 
function does for long times. Based upon the = 2, 4 
examples, the answer to this question might depend upon 
whether N/2 is even or odd. 
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Finally, we turn to the unresolved question, vigorously 
debated in the literature for over a decade, 
of whether the large- N limit of the two-spin correlation 
functions for rings of classical Heisenberg spins, based 
on the dynamics of Eq. (Q) , will decay to zero with the 
leading behavior t~ x / 2 for long times. This asymptotic 
behavior can easily be derived for an infinite linear chain 
of classical spins whose dynamics are governed by the 
following discretized version of a continuous spin hydro- 
dynamics, 



dSi 



•Sj_i-2S,). 



(84) 



The parameter 7 = D/a 2 , where a is the lattice constant, 
and D is a spin diffusion coefficient. The final result 
for the vector component \x of the two-spin correlation 
function for a pair of spins n lattice sites apart coincides 
with the probability distribution for a one-dimensional 
continuous time random walk from the origin, p6fl 



exp(-t*)/„(t*), 



(85) 



where now t* = 2^t and I n (z) is a modified Bessel func- 
tion. For fixed n, the leading behavior for large t is indeed 
proportional to (i*) -1 / 2 . 

For the four-spin ring with diffusive dynamics, Eq. 
(p4|), one finds at infinite temperature, 



and 



but 



C 22 (i) = [l + cxp(-r)] 2 /4, 
Cu(t) = [l-exp(-i*)] 2 /4, 



C 12 (t) = [1 - cxp(-20]/4, 



<S 2 W) = [l + exp(-20] 2 /4 



(86) 
(87) 



(88) 



(89) 



for i = 1, . . . , 4. Although the three correlation functions 
are different, satisfying the conservation law, Eq. (50), 
and C\2{t) correctly reaches its asymptotic limit faster 
than do either of the other two, they all go exponentially 
to the same asymptotic value 1/4, and do not exhibit any 
oscillations. In addition, the rms value of each individual 
spin magnitude decays from its assumed initial value of 
unity to 1/2. 

Generalizing this treatment to a closed finite ring with 
N sites, as derived using Eq. (|84|), leads to the follow- 
ing results. The corresponding time correlation functions 
all approach the same non-zero limit 1/N for long times, 
with correction terms that rise or decay as exp(— n^t*), 
where = 1 — cos(2tt/N). In addition, each rms in- 
dividual spin magnitude decays exponentially from its 
initial value of unity to 1/7V 1 / 2 . These features are dis- 
tinctly different from the exact results we have obtained 
for the four-spin ring, and others have obtained for the 



dimer and three-spin ring, [|25|,[23| based upon Heisenberg 
dynamics, Eq. (Q). However, the diffusive approximation 
does suggest that for N divisible by 4, Ci^+N/^t) should 
approach its asymptotic limit faster than the other cor- 
relation functions, as occurs in the exact treatment for 
N = 4. It would be interesting to see if qualitatively 
similar features are obtained for higher TV values using 
Heisenberg dynamics. 

It is noteworthy that as long as t is sufficiently small 
compared to N 2 /(2'y), the numerical results of the treat- 
ment for the iV-spin ring based upon Eq. ( |34"| ) are vir- 
tually indistinguishable from the result, Eq. (|85|), for 
the infinite chain. This can be expected, since the rms 
distance achieved by the corresponding random walker 
remains small compared to the circumference of the ring, 
hence for all intents and purposes the behavior should be 
the same as that of an infinite linear chain. 

Now our main point is that diffusive spin dynamics 
based upon Eq. ([si]), while conserving the components 
of the total spin vector, does not preserve the length 
of the individual spin vectors. However, this property, 
|Sj(f)| = 1 for i = 1, . . . , N, is of course maintained at all 
times by Heisenberg dynamics, Eq. (^). That is, the spin 
dynamics based upon Eq. $3\ ) are fundamentally differ- 
ent in character from those based upon Eq. Q). Not 
surprisingly, the present exact results at infinite temper- 
ature, Eqs. ( p6| ) - (p8|), displayed in Fig. 3a, along with 
the analytic expansions for long times, Eqs. ( |59| ) - (|6l|), of 
the two-spin correlation functions for N = 4 spins based 
upon Heisenberg dynamics, Eq. (|7|), differ both quanti- 
tatively and qualitatively from the results, Eqs. (^6|) - 
(p9|), that one obtains using diffusive dynamics from the 
four-spin version of Eq. (|84|) . These fundamental differ- 
ences only underscore the fact that the broader question, 
of whether the time correlation functions derived from 
the Heisenberg spin dynamics, Eq. (Q), will show a i -1 / 2 
long-time approach to their asymptotic limits, as does 
Eq. (HH) , is acute and deserving of greater attention. 
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VI. APPENDIX 

Here we list the functions /i(s), gi(s), and hi(s) for the 
correlation functions at infinite temperature. We find, 
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A(s) = —s 2 (s-2) 2 (s + 4) 

(5s 2 + 48s- 150), 



2880 



.9i(*) = %^(5* 3 + 12 5 2 -6 S -8), 



f2(8) 



2880s 
(4-s 2 )r 



12 S + ll S 3 + S 5 ln(^) 
15 {2 + .^ 



480s 2 

-( 12 -^ + ^)i»(|±i) 

/ 3 (s) = ai(s) + 02(s)(2-a)In2 



+ 



s 3 (40-s 2 ) 
480 



In 



is + 2) 



+a 3 (.s)(2-s) 2 ln(4-s 2 ), 



where 



ai(s) = 115 1 20Q (61264 - 28800s + 19160s 2 
-4160s 3 + 45s 4 + 188s 5 ) , 



a2(s) 



and 



1 



(96 + 48s + 608s 2 



1920s 2 V 
+304s 3 - 318s 4 - 23s 5 - 4s 6 ) , 



03(5) = 



1 



320s 2 



(4 + 4s + 13s 2 + 12s 3 



In j 



s 3 (s 2 + 32) 
480 

+b 3 (s) ln(s - 2) + & 4 («) ln(s + 2), 

^(-17280- 28496s + 60000s 2 
+21800s 3 - 5600s 4 - 1485s 5 + 188s 6 

62(s) = -|io§(- 96 + 48s - 608s2 



9 3 {s) = 6i(s)+6 2 (s)ln2 + 
where 

Ms) = 



+304s 3 + 318s 4 - 23s 5 + 4s 6 ) , 
^96 + 96s + 140s 2 



Ms) 

and 

h{s) = ~38W 
and finally, 



(«~2) 2 
3840s 2 

+ 116s 3 + s 4 + 4s 5 



(90) 
(91) 

(92) 



(93) 



(94) 



(95) 



(96) 



(97) 



(98) 



(99) 



(100) 



(s-2)(s + 2) 4 (6-9s + 4s 2 ), (101) 



h 3 (s) = 



(s-2) 
3840 



+2(92 + 46s - 57s 2 ) ln(^-^) 
-10(8 + 4s-3s 2 )ln( 



s + 2 



+ 



48(5 + 2) / (s-2) 2 x 
1 2{s + 2)) 



+ ci(s), 



(102) 



where 



Cl ( s ) = — ^- f 2880 + 16152s - 22668s 2 
v y 230400s V 

+6242s 3 + 267s 4 - 188s 5 ) . 
For the three-spin ring, 



and 



* 2 (5-* 2 ) 
15 



.9 4 (s) = ^r^-(-3-2s + 9s 2 + 4s 3 ) 
120s 



(103) 



(104) 



(105) 



To obtain the leading long-time behavior of the 
lim-r^oo Cij(t), we require the functions fi{s),gi(s), and 
hi(s) and their derivatives at the integration endpoints 
s = 0, 2, 4, 6 (0, 1, 3 for the three-spin ring). The relevant 



(106) 
(107) 

(108) 

(109) 

(110) 

(111) 

(112) 

(113) 
(114) 
(115) 

(116) 

(117) 

(118) 

(119) 

(120) 



quantities are 






/i(0) 


= /i(0) 


= .9i (4) = 


ffl (4) 


/i(2) 


= 5i (2) 


17 

~ 180' 




A (2) 


= 5i(2) 


17 

~ ~360' 




A» 


2 

~ 9' 






A (2) 


= 5i (2) 


43 
~ ~360' 




A"(0) 


3 

~ ~16' 


21 
~ 80' 




A" (2) 


= 9i (2) 




<?i" (4) 


1 

~ ~4' 






MO) 


= A (2) 


= h 3 (6) = 


h' 3 (6) 


A(o) 


1 

~ 6' 






A(2) 
/s(0) 


7 

~30 
3739 

7200 


f — In 2, 
15 

43 , „ 

In 2 

120 


5 


A(2) 


= 53(2) 


7 3 

~ 15 ~ 5 


In 2, 


A(0) 


1 

~ ~4' 






A (2) 


11 
~ 40 ~ 


29, „ 
— In 2, 
30 





17 



' ,^ 21 29 , 

o,(2) = In 2, 

y 40 30 



5s(4) = /i 3 (4) 
and 



69 177 
1600 + "80" 



In 2 



459 
320 



In 3, 



1233 



, , 5 57 

y3V ; 3W 36 20 320 
For the three-spin ring, we require 



/4(0) = / 4 (0)=. O4 (3)=.g 4 (3) = 0, 

/4(1)=S1(1) 1 



In 3. 



5' 



/;(D = - 
.94(1) = - 

and 



15' 
17 

15' 



.94(3) = 1. 
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